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ABSTRACT 

The three-point correlation function of cosmological fluctuations is a sensitive probe of the physics of 
inflation. We calculate the bispectrum. _Bg(fci, fc2, fcs), Fourier transform of the three-point function of 
density peaks (e.g., galaxies), using two different methods: the Matarrese-Lucchin-Bonometto formula 
and the locality of galaxy bias. The bispectrum of peaks is not only sensitive to that of the underlying 
matter density fluctuations, but also to the four-point function. For a physically-motivated, local form 
of primordial non-Gaussianity in the curvature perturbation, $ = ^ -I- /nl0^ + .9nl'/'^, where is a 
Gaussian field, we show that the galaxy bispectrum contains five physically distinct pieces: (i) non- 
linear gravitational evolution, (ii) non-linear galaxy bias, (iii) /nl, (iv) /nl' ^'^^ (^) ffNL- While (i), 
(ii), and a part of (iii) have been derived in the literature, (iv) and (v) are derived in this paper for the 
first time. We also find that (iii) receives an enhancement of a factor of ~ 15 relative to the previous 
calculation for the squeezed triangles (fci « ^2 ^ ^a)- Our finding suggests that the galaxy bispectrum 
is more sensitive to /nl than previously recognized, and is also sensitive to a new term, ^nl- For a 
more general form of local- type non-Gaussianity, the coefficient /^^ can be interpreted as tnl, which 
allows us to test multi-field inflation models using the relation between the three- and four-point 
functions. The usual terms from Gaussian initial conditions, (i) and (ii), have the smallest signals in 
the squeezed configurations, while the others have the largest signals; thus, we can distinguish them 
easily. We cannot interpret the effects of /nl on _Bg(fci, fc2, fcs) as a scale-dependent bias, and thus 
replacing the linear bias in the galaxy bispectrum with the scale-dependent bias known for the power 
spectrum results in an incorrect prediction. As the importance of primordial non-Gaussianity relative 
to the non-linear gravity evolution and galaxy bias increases toward higher redshifts, galaxy surveys 
probing a high-redshift universe are particularly useful for probing the primordial non-Gaussianity. 
Subject headings: cosmology : theory — large-scale structure of universe 



1. INTRODUCTION 

Are primordial fluctuations Gaussian, or non- 
Gaussian? The simplest models of inflation, driven by 
a slowly-rolling single scalar field with the canonical ki- 
netic term originated from the Bunch-Davis vacuum, pre- 
dict the amplitude of primordial non-Gaussianity that 
is below the detectable level. Therefore, a convincing 
detection of primordial non-Gaussianity would rule out 
the above simplest models, and thus lead to a break- 
thro ugh in our understa nding of the physics of inflation 
(see iBartolo eral|[20p . for a review . 

The tightest limit on primordial non-Gaussianity so 
far comes from the angular bispectrum, spherical har- 
monic transf orm of the ang ular three-point correlation 
function (see lKomatsull200ll for a review), of anisotropy 
in t he cosmic microwave background (CMB) radiation 
CseelKomatsu et al.|[200l ISmith et al.]|2009l : ICurto et all 
120091 for the latest limits) . 

The large-scale structure of the universe can also 
provide alternative ways of probing primordial non- 
Gaussianity through abundances and clustering proper- 
ties of galaxies and clusters of galaxies. However, as the 
large-scale structure of the universe is more non-linear 
than CMB, it was generally thought that CMB would 
be the most pro mising way of con straining primordial 
non-Gaussianity (I Verde et al.ll200GD. 

On the other hand. ISefusatti &: Komatsiil (|2007[ ) have 



shown that observations of the large-scale structure 
of the universe in a high-redshift universe, i.e., z > 
1, can provide competitive limits on primordial non- 
Gaussianity, as the other non-linear effects are weaker 
in a high redshift universe. Specifically, they calcu- 
late the bispectrum of the three-dimensional distribution 
of galaxies, B„{k^ , k■2,k:^ )} on large scales as (see also 
IScoccimarro et all (|2004r) ) 



Bg{ki,k2,k3,z) 



2bl ^F^'\h,k2)P,n{ki,z)P„,{k2,z) + (2 cyclic) 
bjb2 [P,n{ki, z)P„,{k2, z) + (2 cyclic)] , (1) 



P,n{ki,z)P,n{k2,z) kjTjks) 

klT{ki) klT{k2) D{z) 



+ (2 cychc) 
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where iJo and Vim are the present-day value of Hubble's 
constant and the matter density parameter, respectively, 
Pm{k, z) is the power spectrum of linear matter density 
fluctuations, D{z) is the linear growth factor, T{k) is the 
linear transfer function whose limit is T{k) — > 1 as fc — > 0, 

and F!f\ki,k2) is a known mathematical function given 

^ The bispectrum, the Fourier transform of the three- 
point correlation function, is defined as (5{fci)(5(fc2)<5(fc3)) = 
(27r)^_B(fci,fc2, fc3)5°(fci + k2 + ks). 
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by (|Bernardeau et al.ll2002f) 

'ih,k2) = - 



fc2 



ki ■ k2 



2fcifc2 \fc2 ki J 7 \ fcife 

' (2) 

This function vanishes in the squeezed limit, ki — ~hi 
(the triangles with the maximum angle, i.e., tt, between 
fci and kfi , and | fci | = | ^2 1 ) , and takes on the maximum 

value, P'^f^ = (a + l)^/(2a) > 2, in the opposite limit, 
k\ ~ ak2 where a > 1 (the triangles with the vanishing 
angle between ki and k2). 

Here, bi and 62 are the linear and non-linear galaxy 
bias parameters, respectively, which relate the underly- 
ing matter dens ity contrast, Sm, to the galaxy density 
contrast, Sg, as (|Frv fc Gaztanagalll993f l 



Sg{x)=bMx) + ^[Sl{x)-a^] + 



(Sf^), which ensures {Sg) — 0. 



(3) 



where a 

The last two terms in Eq. ^ are the well-known results 
for G aussian initial conditions (see iBernardeau et all 
I2OO2I for a review), whereas the first term is the effect 
of the primordial non-Gaussianity of the "local type," 
whose Bardeen's curvature perturbation, is written 
as ^(x) = (l){x) + /NL<^^(=g)i where (j) is a Gaussian field 
jSalopck & Bond" ''1990": 'Gangui et all 1X9941 : IVerde et all 
|2)00: Ko matsu fc S pcrgcl 2001). The latest limit on thi s 
parameter is /nl = 38±21 (68% CL) (jSmith et al.ll2009D . 

However, Sefusatti fc Komatsu's equation, Eq. ([1]), 
may req uire modifi c ations , in light of recent an- 
alytical iD alal et al.' '2008': 'Matarrese fc V erdd l2008l: 
Slosar et al. 2008: Afshordi fc Tolley 2008: Ta ruva. et ail 
2008": 'McDonald 2008) and numerical (Dalai et al .ll2008l: 
Desiacques ct al.,i.200& iPillepich et al.i.2008;,Grossi et al.1 



2OOI studies of the effects of primordial non-Gaussianity 
on the galaxy power spectrum. These studies have 
discovered an unexpected signature of primordial non- 
Gaussianity in the form of a scale- dependent galaxy 
bias, i.e., Pg{k,z) = bl(z)Pra{k, z) [bi(z) + 

A&(fc, z)]^P,„(fc, z), where 



A6(fc,z) 



3(6i(z) - l)/NLa 



iHqSc 



D{z)k^T{k) 



(4) 



and Sc ~ 1.68 is the threshold linear density contrast for 
a spherical collapse of an overdensity region. 

Then, several questions arise: can we still use Eq. ([T]) 
for the bispectrum? Should we replace bi by 61 -I- A6(/c)? 
Does the first line in Eq. ([1]) somehow give the same 
correction as Ab{k)7 How about 62'? We are going to 
address these questions in this paper. 

2. BISPECTRUM OF DARK MATTER HALOS 

In this section, we derive the galaxy bispectrum for 
non-Gaussian initial conditions by using two different 
methods. In ij2.11 we shall use the "functional integra- 
tion method" for computing n-point correlation func- 
tions of peaks of the cosmological density fluctuations 
(jPolitzer fc Wise 1984: Grinst ein fc Wise 1986). In g221 
we shall present an alternative derivation of the same 
result by using a local bias assumption. 

2.1. Mararrese-Lucchin-Bonometto (MLB) method 



We sh all use the Matarrese-Lu cchin-Bonomctto (MLB) 
formula (jMatarrese et al.lll986[ ) which allows one to cal- 
culate the n-point correlation functions of peaks for 
non-Gaussian initial conditions. Th is approach is espe- 
cially well suited for our purposes, as lMatarrese fc Verdd 
(|2008f l have applied the MLB formula to compute the 
scale-dependent bias of the galaxy power spectrum. We 
shall apply the MLB formula to compute the galaxy bis- 
pectrum for general non-Gaussian intial condition. 

We study the three-point correlation function of the 
spatial distribution of dark matter halos. Let us con- 
sider the probability of finding three halos within three 
arbitrary volume elements: dVi, d Vp, , and dV^ , w hich are 
at Xi, X2, and x^, respectively, as (|Peebleslll980() 

P{xi,X2,X3)^n^ [1 + Oii^u) + S.h{x23) + S.h{x3l) 

+ Ch{xi,X2,X3)]dVidV2dV3, (5) 

where Xij = \xi — Xj\, and and are the two- and 
three-point correlation functions of halos, respectively. 

The next step is to relate the correlation functions of 
halos, and (h, to those of the underlying matter dis- 
tribution function. The locations of halos coincide with 
those of the peaks of the matter density fluctuations; 
thus, one can compute S,h and Ch by computing the corre- 
lation functions of peaks above a certain threshold, above 
which the peaks collapse into halos. 

We shall assume that halos would be formed in the re- 
gion where the smoothed linear density contrast exceeds 
Sc- For a spherical collapse in an Einstein-de Sitter uni- 
verse Sc = 3(127r)^/'^/20 ~ 1.68, and one can find other 
values in the ellipso idal collapse in arbitra ry cosmologi- 
cal models (see, e.g.. lCoorav fc Shethl[200l . for a review). 
The mass of halos is determined by the smoothing radius, 
R, i.e., M = (47r/3)pm-R'^, where pm is the average mass 
density of the universe. The smoothed density contrast, 
is related to the underlying mass fluctuations, 5m, 
as Sii{x) = J d'^a/WB.{\x ~ a/|)(5m(a/), where Wr{x) is a 
smoothing function. We shall use a top-hat filter with 
radius R for Wii{x). 

Using the MLB formula, we find 

1+^^(2:12) + ^hix23) + S.h{x3l) + Ch{xi, X2, X3) 



-exjp 



(2), 



-3 l,mi— 1712—0 



Xi, 



, Xi,0C2, ■ 



, X2, X3, 



, X3 



mi times 7712 times 7713 times 



-3 



^1 



77 times 



(6) 



where 7773 = 77—7771 ^'^2, 1^ = Sd (Jr, <j% is the variance of 
matter density fluctuations smoothed by a top-hat filter 

(n) 

with radius R, and denotes the connected parts of 
the 77-point correlation functions of the underlying mat- 
ter density fields smoothed by a top-hat filter of radius 
R. Here, we have assumed that we are dealing with high 
density peaks, i.e., ly ^ 1, which arc equivalent to highly 
biased galaxies, 61 ^ 1. 

As ^j^'' ^ 1 on the large scales that we are interested 
in, we expand the exponential in Eq. We keep the 
terms up to the four-point function, as this term provides 
the dominant contribution to the three-point function. 
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Wc find 

Ch{xi, X2, X3) = -^^^^^(a;i, X2, X3 



+ 



R. 



24 



eR'{^i2)CR'{x2z) + {2 cyclic) 



£.r\xi-, Xi,X2, Xz) + (2 cyclic) (.7) 



The bispectrum of halos in Lagrangian space, 
Bl^{ki, ^2,^3), is the Fourier transform oi (h{xi, X2, x^): 



B^(fcl,fc2,fc3) 

;/3 r 

= -3- Bn{ki,k2,k3) 



'R L 



— {PR{ki)PR{k2) + (2 cyclic)} 

(^R 



Tniq, ki - q, fe, fe) + (2 cyclic) 



(8) 



where Tr is the trispectrum, Fourier transform of ^j^^ . 
Here, we call Bf^ the Lagrangian space bispectrum, as it 
relates the halo over-density to the initial matter over- 
density with its amplitude extrapolated to the present 
epoch. If we assume that the halos move in the same way 
as matter, the observed bispectrum in Eulerian space, 
Bh, would be the same expression with Eq. ([8]), except 
for the coefficients: 



= bl 



Bh{ki, k2, ks) 
BR{ki,k2, h) 



bi 



{PR{ki)PR{k2) + (2 cyclic)} 



24 



(27r)' 



7Tii{q, ki - q, k2, k^) + (2 cyclic) 



(9) 



Here, hi is the so-called linear Eulerian bias parameter, 
5i = 1 + vjaj^^ and &2 is the non-linear bias parameter. 

2.2. Alternative derivation 

In this section, we present an alternative derivation of 
the galaxy bispectrum, Eq. ([5]). On large enough scales, 
wc may approximate the relation between the galaxy dis- 
tribution and the underlying density fiuctuation as a lo- 
cal function. We then Taylor-expand this local function 
in a power series of (see Eq. ([3])). 

When computing the correlation functions of halos of 
a given mass M, we may smooth the matter density field 
with the same filter over the corresponding length scale 
R, Wr{\x— a/ 1), which was defined in the previous sec- 
tion. We then Taylor-expand Eg in a power series of the 
smoothed density field, Sr{x)^ as 



6g{x) = h^bR{x) + ^ [4(a;) - • • • 

In Fourier space, one finds 
5g(k) = hr8R(k) 



(10) 



&2 

" 2 



q 



i2ny 



5R{k^q)5R{q)~cj'J''{k) 



(11) 



^ Note that this expression, 61 = 1 + v/ajj., ag r ees w ith the 
Hnear halo bias parameter derived by IMo &: White! 1)19961 ) , 61 = 
1 + {u'^ — l)/5c, for high density peaks, fc" 3> 1. 



where 5^ is the Dirac delta function. Wc calculate the 
bispectrum of galaxies directly from Eq. pip : 

{5M5g{k2)5g{h)) 
^bl{SRih)dR{k2)6Rik3)) 



bib2 



(27r)= 



(<5fl(fci - q)SRiq)SRik2)6R{k3)) 



<jj,6''ih){5R{k2)dRih)) + (2 cyclic) 



(12) 



The first term of Eq. (1121) is the matter bispectrum, 

{SR{h)6R{k2)SR{k3)) - (27r)3BflXfci,A:2,fc3)<5^(fci23), 

where fci23 = fci + fc2 + ^3. Wc further calculate the en- 
semble average of the four-point function in the second 
term of Eq. (|12p . For non-Gaussian density fields, four- 
point function is given by a sum of connected (trispec- 
trum) and unconnected (products of the power spectra) 
parts as 

{SRih - q)SR{q}SRik2)6R{k3)) 
= (2^)6Pfl(q)Pfl(A:2)(5^(fci)5^(fc2 + h) 

+ (27r)6p«(fc2)FK(fc3)5^(fc2 + q)S''ih + h - q) 
+ (27r)6pfl(fc2)Pfl,(fc3)'5^(fc3 + q)S''{k2 + ki - q) 
+ (27r)3Tfl(fci - q, q, fca, fc3)'5^(fci23), 

where Tr is the matter trispectrum. Note that the 
first term in the above equation cancels the last term 
in Eq. (fT2|l . Combining the above equations, Eq. p2|l 
becomes 



Bh{ki, ka, ks) 
(2^) 



Pfi(fci, fc2, fc3) + 7^ {PRiki)PRik2) + (2 cyclic)} 



1 b 

'2b 



- I 7?T?^«(9, fci - fc2, ^3) + (2 cyclic) 

1 J [21^ r 



.(13) 



Therefore, we find that the MLB method and the local- 
ity bias assumption give formally the same results. This 
may imply that the process of picking up density peaks 
over a certain threshold is the local process. 

At the same time, there is a subtle difference between 
Eq. ([9]) and Eq. p3)) : the coefficient in the last term. 
By evaluating the last cyclic term in Eq. for the lo- 
cal type of non-Gaussianity, we find that the integration 
of the smoothed trispectrum depends on the smoothing 
scale, i?, up to a constant factor of l/<j\ on large scales, 
say, k < 0.1 h/Mpc. For example, the bottom right panel 
of Fig. [n] shows that B^^ and B^^ , which are defined 

in Eq. (|29)l and Eq. ((30)l . respectively, do not depend 
on the smoothing scale R, as they include I/ctr in their 
definitions. 

Therefore, it is physically more sensible to include 
explicitly in the equation such that the dependence on 
smoothing scales on large scales can be absorbed by the 
bias parameters. This motivates our writing the final 
form of the halo bispectrum, derived from the local bias 
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assumption, as 
Bh{ki, fe, fca) 



BflXfci, fc2, h) + ^ {Pfl(fci)Pfl(fc2) + (2 cyclic)} 



'-r^TT / 77r^TR{q, h - q, fe, fe) + (2 cyclic) 
Ol 2cr^ J [Ztt)-^ 



(14) 



with three bias parameters: &i, &2 and 62 = ^2f/;- Note 
that 62/61 — * (5c for the MLB formula. Although Sc is 
known to be 1.68 for the spherically collapsed halo in the 
flat matter dominated universe, its precise value, in this 
context, needs to be tested against N-body simulations. 

Eqs. Q and (fT4|) are the first main results of this pa- 
per, which are general and can be applied to any mod- 
els of non-Gaussianities, once we know the bispectrum 
and trispectrum of the underlying matter density field. 
Note tha t Eg. ([T4|) was also obtained independently by 
iSefusattil ((20091) . 

In principle both 61 and 62 are calculable from the the- 
ory of collapse of halos (see, e.g., ICoorav fc Shethir2002l . 
for a review); however, in practice it is more convenient 
and safe to treat them as free parameters that one should 
marginalize over when ext racting the cosm ological infor- 
mation, such as /nl- See I Jeong fc Komatsu (2009. ) for 
the same argument applied to the galaxy power spec- 
trum. 

3. EFFECTS OF LOCAL-TYPE PRIMORDIAL 
NON-GAUSSIANITY ON THE HALO 
BISPECTRUM 



In this section we shall evaluate Eq. (|T4|) for the local- 
type primordial non-Gaussianity with a high-order term 
added: 



$(a;) 



c) + /nl [^^(x) - (02)] + gNL^Hx). (15) 



The cubic-order term does not generate the bispectrum 
of CMB anisotropy or the matter density fluctuations at 
the leading order; however, it does generate the trispec- 
trum, and the CMB t r ispectrum has been calcu lated by 
lOkamoto fc Hul(|2002r ): lKogo fc Komatsul (|2006[ ). On the 
other hand, the bispectrum of halos receives the contribu- 
tion from the trispectrum (see the last term in Eq. (|14p ). 
and thus it is necessary to include 5nl- 

To calculate various components of the bispectra 
shown in Eq. (jl4p . we calculate the transfer func- 
tion, T{k), and the power spectra with the cosmolog- 
ical parameter s in T able 1 ("WMAP-^BAO-fSN") of 
iKomatsu et"al] (|2009f) . 

As for the smoothing scale, we use R = 1 /i^^ Mpc. 
Although the smoothing scale explicitly appears in the 
equation, it makes negligible differences for the bispec- 
trum on large scales, k <^1/R. 

Note that we shall adopt the non-standard convention 
in which $(a;) is Bardeen's curvature perturbation ex- 
trapolated to the present epoch, z = 0, using the lin- 
ear growth factor of <&, g{z) = (1 -I- z)D{z). Therefore, 
our /nl and ^nl in this paper are different from those 
in the literature on the CMB non-Gaussianity by a fac- 
tor of g(1090)/5(0), i.e., /nl = [ff(1090)/5(0)]/CMB and 



5NL = [5'(f090)/.92(0)]5gMB 3, 
The bispectrum and trispectrum of $ are given by 



B^ih, h2, h) = 2/nl [^0(^1)^0(^2) 
and 



(2 cyclic)] , (16) 



r$(fci, ^2, ks, ki) 

= 65NL [P4ki)P<p{k2)P^{h) + (3 cyclic)] + 2/^^ 

X [P^{ki)P^ik2) {P^kis) + P0(fci4)} + (11 cyclic0|7) 

with kij ~ \ki + kj\, respectively. 

While Eq. (fTT]) is the consequence of Eq. (fT5|) , general 
multi-field inflation models do not necessarily relate the 
coefficients of the trispectrum to that of the bispectrum. 
Therefore, one may generalize Eq. (|17p by replacing /^^ 
with a new parameter, tnl, which may or may not be 
related to /nl: 



r$(fci, A;2, fcs, k4 



25 



= 65nl [P4,{ki)P^{k2)P^{k3) + (3 cyclic)] + — tnl 

lo 

X [P4,{k,)P^{k2) {P^ikis) + P4,{ku)} + (11 cycli 

Note that the coeffici ent of tnl reflects t he de finition 
of TNL introduced by iBoubekeur fc LvthI l|2006( l. This 
opens up an exciting possibility that the galaxy bispec- 
trum can test whether tnl = (6/nl/5)2 or other predic- 
tions for the relation between tnl and /nl are satisfied 
observationally. 

We transform these spectra to those of the smoothed 
linear density contrasts, using the Poisson equation. 



2k^T{k) ~ 



3 -//q 



VFfl(fc)$(fc) = A^i^(fc)$(fc), (19) 



where Wi?,(fc) is the Fourier transform of the top- hat filter 
with radius R. Note that Wii{k) — > 1 as fc ^ 0. In 
general W^/f^fc) ~ 1 for fc <C 1/i?. Then, we calculate the 
n-point function of the matter density fields from the 
corresponding correlator of curvature perturbations by 

n 

(4^)(fci) • • • 4'^(fc„)) = Y{MB{h){Hki) ■ ■ ■ $(fc,o). 

1=1 

3.1. Known Terms 
3.1.1. Formula 

The first term in Eq. ([T^ contains the bispectrum 
of matter density fluctuations, which consists of 

two pieces: (i) the non-linear evolution of gravitational 
clustering (B^) and (ii) primordial non-Gaussianity 

(5„G0)4. 

^ /NL ' 

Bniki, k2,k3) = Bg{ki, fc-2, fca) + hi^B]^^^{ki, fc2, fca), 

(20) 

^ The ratio 3(1090)/g(0) is 1.308 for t he cosmological parame- 
ters in Table 1 ("WMAP-I-BAO+SN") of lKomatsu et al.l ||2009|) . 
We ignore the following term in Bfl(fci, k2,k3): 

n "^Rih) I -^F^^'\q, fci - q)T{q, fci - g, fca, ^3) + (2 cyclic), 

where T is the unfiltered primordial trispectr um which contains 
/^j and (/NL- This term is negligibly small IjScoccimarro et al.l 
[20M) . 
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where 

fc2, fcs 



iWR{ki)WB.{k2)Wnik3)2F^'\ki,k2) 
xPm{ki)Pm{k2) + (2 cyclic), (21) 



with f!^'' given by Eq. ([2]), and 

B]^^{k,M,h) 

3 

= 2\{MR{h) [P4ki)P^{k2) + (2 cyclic)] 

i=l 

= 2#^#%\-^«(^3) + (2 cyclic). (22) 
MB\ki) MR(k2) 



One finds that Eqs. ([2^ and (PT|) agree with the first and 
the second terms in Eq. ([T]) on the scales much larger than 
the smoothing scale, i.e., k <^ for which Wr 1. 

One might wonder if it is OK to include the bispcc- 
trum from non-linear evolution of density fluctuations in 
the MLB formula, as Eq. ^ is usually used for the La- 
grangian density fluctuations, i.e., "initial" fluctuations. 
However, it is perfectly OK to use the evolved density 
fluctuations in this formula, as one can always use the 
evolved density fluctuations as the initial data. For ex- 
ample, we can think of starting our calculation at z = 10, 
and ask the MLB formula to take the initial condition at 
z = 10, including non-linear correction. Since wc know 
how to compute the bispcctrum, trispectrum, etc., of the 
underlying mass distribution at z = 10 (including non- 
linear effects), we can use this information in the MLB 
formula. In other words, the "initial" distribution does 
not need to be primordial. We can provide the evolved 
density field as the initial data, and compute the peak 
statistics. The MLB formula does not care whether the 
source of non-Gaussianity is truly primordial or not: the 
only conditions that we must respect for Eq. ^ to be 
valid are (i) high peaks (i/ 3> 1), and (ii) n-point correla- 

(n) 

tion functions are much less than unity, 1, so that 

the exponential in Eq. ([6]) can be Taylor-expanded. In 
this case one would lose an ability to calculate the bias 
parameters, bi and 62, using, e.g., a halo model; how- 
ever, this is not a disadvantage, as the halo model cal- 
culations of the galaxy bias parameter, bi and 62, are at 
best qualitative even for Gau ssian initial conditions (see, 
e.g.. lJeong fc Komatsul[2009[ l. In our approach the coef- 
ficients of individual terms in Eqs. ^ and p4|) . including 
Sc, are free parameters, and need to be determined from 
observations themselves. 

3.1.2. Shape Dependence: Non-linear Gravitational 
Evolution and Non-hnear Galaxy Bias 

How about the shape dependence? First, let us re- 
view the structure of B^{ki, k2, k^) (Eq. pT|) ). which 
has been stud ied extensively in the literature (see 
iBernardeau et a l. 2002, for a review). 

Here, let us offer a useful way of visualizing the shape 
dependence of the bispcctrum. Wc can study the struc- 
ture of the bispectrum by plotting the magnitude of B^^ 
as a function of ^2/^1 and k^/ki for a given fci, with 
a condition that fci > ^2 > k^ is satisfied. In order 
to classify various shapes of the triangles, let us use 
the following names: squeezed (fci ~ fc2 ^ ^3), elon- 
gated (fci = fc2 + fca), folded (fci = 2fc2 = 2fc3), isosceles 



(fc2 = fca), and equilateral (fci ~ k2 = k^). See (a)-(e) of 
Fig. [1] for the visual representations of these triangles. 

The top-left panel of Fig. [5] shows B^ for fci = 
0.01 h Mpc~^. In this regime PR(fci) is an increas- 
ing function of fci (recall that Pnik) peaks at fc w 
0.02 h Mpc~^). Let us then pick the first term in 

Eq. (I21|), F^'\ki, k2)PR{ki)PR(k2), and ignore the cyclic 
terms for the moment. It follows from Eq. ^ and the 

descriptions below it that ^2^*' (fci, ^2) vanishes in the 
squeezed limit (fci = — fc2) and reaches the maximum 
in the opposite limit (fci = ak2). Therefore, we would 
expect this term to give large signals in the "elongated 
configurations," fci ~ k2 + k^; however, as PR(fc) at 
fc < 0.02 h Mpc~^ is an increasing function of fc, one 
can also get large signals when fci and fc2 are equally 
large, fci = fc2. As we have zero signal in the squeezed 
limit, fca = 0, it follows that we can find a large signal 
in the equilateral configuration, fci = fc2 = fca. A similar 
argument also applies to the cyclic terms. As a result, 
for fci = 0.01 h Mpc^^, we find the largest signal in the 
equilateral configuration, and then the signal decreases 
as we approach the squeezed configuration, i.e., the sig- 
nal decreases as we go from (e) to (a) in Fig. [1] 

The top-right panel of Fig. [2] shows B^ for fci ~ 
0.05 h Mpc"\ In this regime PR{ki) is a decreasing 
function of fci, and thus the equilateral configurations 
arc no longer as important as the folded ones. Instead 
we have large signals in the folded configurations as well 
as in the elongated configurations. Note that the exact 

(s) 

squeezed limit is still suppressed due to the form of F2 ■ 
In summary, B^ usually has the largest signal in the 
folded and elongated (or equilateral, depending on the 
wavenumber) configurations, with the squeezed configu- 
rations suppressed relative to the others. The suppres- 
sion of the squeezed configurations is a generic signature 
of the causal mechanism such as the non-linear gravita- 

(s) 

tional evolution that F2 describes. 

The bispectrum from the non-linear bias term, the sec- 
ond term in Eq. ([9|), has the same structure as B^, but 

it docs not contain i^2*''. As a result the non-linear bias 
term does not have as much suppression as B^ has in the 

(s) 

squeezed configuration. In addition, as F2 enhances the 
signal in the elongated configurations, the non- linear bias 
term does not have as much enhancement as B^^ has in 
the elongated configurations. These properties explain 
the bottom panels of Fig. [2l 

As B^ and the non-linear bias term contain products 
of PR{ki)PR(k2) and the cyclic terms, it is oft en more 
convenient to deal with Q/i(fci, fc2, fc3) given by (jPeebled 
[T980h 

n 11 1 I. \ - Bhiki,k2,k3) 

Qh{ki,k2,k3) = n , p , .„ TT-r, (23) 

PR{ki)PR{k2) + (2 cyclic) 

to reduce the dependence on the shape of the power spec- 
trum. This combination is constant and equal to 6162 for 
the non-linear bias term (see the second term in Eq. ([9])). 

The left and right panels of Fig. [3] show 
BG(fci,fc2,fc3)/[Pfi(fci)Pfi(fc2) + (2 cyclic)] for 
fci = 0.01 h Mpc"^ and 0.05 h Mpc~\ respectively. 
We find that Qh better refiects the shape dependence 

(s) 

of F2 irrespective of fci: it has the largest signal in 
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the folded and cfongated configurations in both large 
and small scales. The squeezed configurations are still 
heavily suppressed relative to the others. 



with r^*^' given by 



3.1.3. Shape Dependence: /nl Term 

How about the /nl term, S]f^^°(fci, fc2, fe)? This 
term has a completely different structure. Let us pick 
the first term, MR{ki)P^(ki)MR{k2)P4,{k2)MR{k'i), in 
Eq. (P^ . The important point is that the power spec- 
trum of (j) is always a decreasing function of k, i.e., 
Pcf,{k) oc 1/k^ for a scale-invariant spectrum. On the 
other hand, on large scales we have T{k) — > 1 and 
7Wi?(fc) oc fc^. Therefore, collecting all the cyclic terms, 
we find B}?f^"(fci,fe,fc3) cx fc|/(fcifc2) + kl/ik.k^) + 
^1/(^2^3) = {kf + k2 + kl)l{kik2k^). In other words, 
it has the largest signal when one of fc's is very small, 
i.e., the squeezed configurations, which is opposite to the 
structures of and the non-linear bias term. The mid- 
dle panels of Fig. H show B^j^^ for ki = 0.01 h Mpc"^ 

and 0.05 h Mpc~^, and we find the largest signals in the 
squeezed configurations. We also find that Qh from the 
/nl term, B]?GO(fc^^ fc^, fc3)/[-Pi?(fci)-Pi?,(fc2) + (2 cyclic)], 
still has the largest signal in the squeezed configurations. 

These properties allow us to distinguish between the 
primordial non-Gaussianity and the other effects such 
as th e non-linear gravitational ey olution and non-linear 
bias. ISefusatti fc Komatsul ()2007f ) have studied in detail 
how well one can separate these effects using Qh- 

3.2. New Term 
3.2.1. Formula 

Now, we shall evaluate the new term, the third 
term in Eq. (|14p. which was not considered in 
ISefusatti fc Komatsul ([2007). The trispectrum is gen- 
erated by the primordial non-Gaussianity, as well as 
by the non-linear evolution of matter density fluctua- 
tions. The non-linear evolution of matter density fluctu- 
ations on large scales is given by perturbation theory (see 
iBernardeau et ai1l2002l for a review) . Let us expand the 
filtered non-linear matter density field in Fourier space 
as 

5R{k)^WR{k) \5^^\k) + 5^^\k) + 5^'^\k) + ---] , 

(24) 

where J^"^(fc) is the n-th order quantity of the linear 
density contrast, 5^^\k). Then, the connected matter 
density trispectrum is given by 

TR{ki,k2, fc3, ^4) 
= T}l^\ki,k2, fc3, ki) + {ri^"'(fci, fc2, fcg, fc4) + (3 cyclic)} 
+ {Ti"3(fci, fc2, A;3, ki) + (3 cyclic)} 
-t- {T],'^\h,k2, h, ki) + (5 cyclic)} + 0{cj)^), (25) 



(2^)^<5^(fcl234)T;^"(fcl,fc2,fc3,fc4) 

4 

n WR{kr,){5^'^ {ki)5^^^ (fc2)(5W (fc3)<5(') {ki)),. (26) 



The leading contributions of all the terms shown in 
Eq. (pS)) arc order of 0^. 

The first term, T^^^, is the linearly evolved primor- 
dial trispectrum calculated from Eq. p7|) . and thus it 
contains the terms proportional to /^^ ^^'^ ffNL- The 
second term, T^^^^, has a coupling between the primor- 
dial non-Gaussianity (linear in /nl) and the non-linear 

gravitational evolution (linear in ^2^"*)- These two terms 
are important on large scales. 

The other terms, T^^^ and T^^^, do not have contri- 
butions from /nl or (/nl at the leading-order level, but 
solely come from the non-linear gravitational coupling; 
thus, they may be ig nored on large s cales we are con- 
sidering in this paper. ISefusatti! (|2009l ) also derived and 
studied T^^^ as well as r^"^ ^nd r^i22^ 

Therefore, we approximate the integration in the third 
term of Eq. as 



2^x1 



24 



[TRiq, h - q, fc2, ^3) + (2 cychc) 



(27r)3 



r]j'^(g,fci-9,fc2,A;3) + (2 cyclic) 



(27r)3 



T^Riq, fci - g, fe, h) + (2 cychc) \ (27) 



where "cyclic" denotes the cyclic combinations 
of fei, fc2, and ks, and and t1^' denote 



Tll\h,k2,k3,ki) 

T^^\ki,k2,k3,ki) 



T^"i(fci,fc2,fc3,fc4), and 
T}l^^{ki,ka,h,ki) + (3 cyclic). 



- R^ y^-Li "-zi "-ii "-i) — ^R 
respectively. 

The first term in Eq. (P7)) is the integration of the lin- 
early evolved primordial curvature trispectrum, which 
contains two pieces: one proportional to /^^ and an- 
other to (7nl (see Eq. ([T7|) ). Therefore, we symbolically 
write the first line in Eq. ((27)) as 



24 



(2^ 



Tll\q,ki- q,k2,h) + i2 cychc) = gNL^"": (^i , , fc; 



fLB]^Akl,k2,k3), 



(28) 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



k /k k /k 

Fig. 2. — Shape of the bispectrum, B(ki ,^21^3)- Each panel shows the amplitude of the bispectrum as a function of k2/ki and fcs/fei for a 
given ki , with a condition that kg < k2 < ki is satisfied. The amplitude is normalized such that it is unity at the point where the bispectrum 
takes on the maximum value. For the visual representations of the triangle names such as the squeezed, elongated, folded, isosceles, and 
equilateral, see Fig. [T] (Top Left) The bispectrum from the non-linear gravitational evolution, (Eq. I I21I I). for fei = 0.01 h Mpc~^. 
(Top Right) for ki = 0.05 h Mpc~^. (Bottom Left) The bispectrum from the non-linear galaxy biasing, PR(ki)Pji[k2) + (2 cyclic) 
(the second term in Eq. (O), for fci = 0.01 h Mpc~^. (Bottom Right) PR{ki)Pii(k2) + (2 cychc) for fci = 0.05 h Mpc~^ 
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.squeezed 
folded 




Fig. 4. 



0,4 0.6 

Same as Fig. [2] but for the terms proportional to /nl 



0.4 0,6 

(Top) the BJ^f term (Eq. (|33]l), (Middle) the 3"^^° term (Eq. iTIl). 



and (Bottom) the BJ^^ term (Eq. (I39l l). Note that the non-Gaussian terms diverge in the exact squeezed limit, ^3 — > 0; thus, we show 

these terms normalized to be unity at fca/fci = 10~^. In order to facilitate the comparison bettor, we draw the dotted contour for all six 
panels. 
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-^quiloteroy 



1.0 
0.9 
^' 0.8 
^"0.7 

0.6 
0.5 









k^ = 0.05[h/Mpc] i 




I— 1 1 



10'' 
10-^ 
10'^ 
10"^ 
10'^ 



Fig. 5. — Same as Fig.|2] but for (Top) the ijnl term (Eq. i30l ). and (Bottom) the /^j^ term (Eq. JSTJ). Note that the non-Gaussian 
terms diverge in the exact squeezed Umit, fca — » 0; thus, we show these terms normaUzed to be unity at fcs/fci = 10^'^. In order to facilitate 
the comparison better, we draw the dotted contour for top panels. 

where 



JNL /n't. 



^MR{k2)MR{k:i) 



(27r)= 



MB{q)MR{\ki - q\)P^{q) 



2<jj 

X [P^{\ki - q\)P^{\k2 + q\) + P^{\ki - q\)P^{\h + q\)] + (2 cyclic) 



'SMR{k2)MR{k^)P^{k2] 



:MR{q)MR{\ki - q\)P^{q)P^{\k:i + q|) + (2 cyclic) 



+ ^MR{k2)MR{k^)P^{k^) J ■^MR{q)MR{\h - q\)P4q)P4\h + q|) + (2 cyclic) 
+ 8MRik2)MRik3)P4ki) [P^ik2) + P^ik^)] [ -^MR{q)MR{\ki - q\)P4,{q) + (2 cyclic) 



'^MR{k2)MR{k^)P^{k2)P^{k:i) 



(2^ 



(2^)3 



MR{q)MR{\ki - q\) 



B^^^{kuk2M) 



X [P4>i\k2 + q\) + P4,{\k3 + q\)] + (2 cyclic 
1 



24 



GMR{k2)MR{k3) [P4k2) + P4k3)] 



d^q 

(2^ 



(29) 



MRiq)MR{\h - q\)P^{q)P^{\ki - q\) + (2 cyclic) 



l2MR{k2)MR{k3)P^{k2)P4,{k:i) j .f3^MR{q)MR{\k^ - q\)P^{q) + (2 cyclic) 



(30) 



I 

We find that the first three cyclic terms in Eq. ([29|) arc for k < 0.1 h Mpc""'^. Therefore, one may just calculate 
parametrically small on large scales and may be ignored the last two cyclic terms: 



1 

24 



8MR{k2)MR{k3)P^{k^) [P^{k2) + P0(fc3)] 



d3q 



-^MR{k2)MR{k3)P^{k2)P^{k:i 

X [P^(|fc2 + q|)+P0(|fc3 + g|)] 



(27r)3 
(2 cyclic 



(27r)3 

MR{q)MR{\k^ - q\) 



MR{q)MR{\ki - q\)P^{q) + (2 cyclic) 



(31) 
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Next, the second term of Eq. (|27|) contains a cross- 
correlation between the non-hnearly evolved density field 
(,5(2) -fi''^ [i^^^-*]^) and the primordial bispectrum, and 

(s) 

thus it is linearly proportional to /nl and . We 



present the explicit functional form of Tj^^^^ as well as 
the full expression of the second term of Eq. ([27]) in Ap- 
pendix]^ Here, we only show the final result. We write 
it as 



1 f 



/NL 



B^{k^,k^, h) + B]^^\h,k2, fcs) + ^B]^^>ik,,k2, fcs) {^^(fci) + 5B{k2) + GRih)} 



(32) 



where 



Bl^{kuk2, k3)=4WRiki)WR{k2)WR{k3) 



24, 



8T^fl(fc2)W^fl(fc3)A^(fc3)-Pm(fc2) 



Pm{ki)P,nik2)F^'\h,k2) + (2 cyclic) !>, 



- q\)WRiq)M{\h - q\)M{\k2 + q\)F^'\-k2, fcs + q) 



X {P0(fc3)P0(|fei - q\) + P4,{k3)P^{\k2 + q\) + P^{\h - q\)Pd\h2 + q\)} + (5 permutation) 



~8WR{k2)WRik3)Mik3) 



1 tJ 



WR{\ki - q\)WR{q)M{\ki - q\)M{\k2 + q\)Pra{q)F!^"\- q. k2 + q) 



(33) 



{P^dfci - q\)P4,{k3) + P^{\ki - q\)P4\k2 + q\) + P^[kj,) P^{\k2 + q\)} + (5 permutation) 



+ 8{WR{k2)WR{kz)] P^{k3)M{k2 



(27r)3 



Mip)Mi\k2 ^ pDPM 



{P^i\k2 - p\) + 2P^{k2)} F!f\p, k2-p) + (5 permutation) 



(34) 



Here, M{k) = MR{k)/WR{k). 



r 



In the above equations we have defined two functions, 
J-R{k) and QR{k), which arc given by 



•^^W^^/ -^PMMR{q)MR{\k-q\\ 



P^{\k-q\) 

pm 

d\ WR{q)WR{\k- 



2-1. 



(27r)3 
xP^{q)F^'\k, 



WR{k) 



(35) 



(36) 



As shown in Fig. [S] and Fig. El both TR{k) and ^flXfc) 
are almost constant on large scales. If we do not have a 
smoothing, i.e., i? — > 0, the large scale asymptotic value 
oiGR{k) is 17/42. However, the presence of filter changes 
this asymptotic value. As k 0, 



GRik) 



13 

84 



44 



WRiq)P^iqf^, (37) 



(27r) 



qR 



whose value depends on the smoothing scale, R, as shown 
in Fig. il 



Let us study the structure of each term in Eq. ([32|l . 
The first piece is B^^ . On very large scales, where 
WR{k) — > 1 and TR{k) 1, B^ becomes a product of 
the usual matter bispectrum for Gaussian initial condi- 
tions, i?^, and the scale dependent bias shown in Eq. 
as 

2/NL-^fl(fc) _ 3/NL^^o^m pR{k) ifni^H^i^m 



MR{k) 



k^T{k) WRik) k^T{k) 



(38) 

as fc — > 0. Therefore, we can interpret this term as a scale 
dependent bias multiplying the usual matter bispectrum 
for Gaussian initial conditions; however, this behavior is 
not generic - in fact, the other terms cannot be expressed 
in terms of products of the scale-dependent bias and the 
results in the continuous limit, Eq. ([T]). 

The next piece is B^^^{ki, k2, fcs). By numerically cal- 
culating Eq. (|34p . we find that the terms that contain 

F2''\q,k ~ q) are parametrically small on large scales, 
and that the dominant contributions come from the first 
permutation terms. Therefore, we approximate Eq. (j34p 
on large scale {k < 0.1/iMpc~^) as 
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24 



[P<pi\k2 + q\) + P<t.i\ki - q\)] F.^''\-k2, k2 + q) + (5 permutation) 



'^''^ Wnilh - q\)WR{q)M{\h - q\)Mi\k2 + q\) 



(2^) 



(39) 





: R 


= 1 Mpc/h / / / / - 
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= 2 Mpc/h ; / / / 




: R 


= 5 Mpc/h / 1 : 1 




: R 


= 10 Mpc/h / / / / 

■' ' ■ / 
' ' ■ / 











0.0001 0.0010 0.0100 0.1000 1.0000 

wovenumber k [h/Mpc] 

Fig. 6. — Shape of the function, J^fi{k), defined in Eq. l|35|l . 
We show Tii(k) for four different smoothing lengths: i? = 1, 2, 5, 
10 Mpc/h. 
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wovenumber k [h/Mpc] 



Fig. 7. — Shape of the function, Gnik), defined in Eq. I I35I I. 
We show G[i{k) for four different smoothing lengths: ii = 1, 2, 5, 
10 Mpc/h. 

How about the last term of Eq. ([5^ . 

4ByG0(fc^, fc^, h) {Giiiki) + gR{k2) + QR{k^)}i As 

Gnik) — > constant on large scales (Fig. [7]), this piece 
becomes B^^^ multiplied by a constant factor whose 
exact value depends on the smoothing scale, R (Fig. [8]). 



In summary, we have derived the new terms in the 
galaxy bispcctrum, which arise from the integration of 
the matter trispcctrum. While we find one term, B^^, 
includes the scale-dependent bias which appears on the 
galaxy power spectrum, we also find that there are more 
terms contributing to the galaxy bispcctrum. 

Eq. ^ along with Eqs. and Eq. ^ along 

0.40 E — ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' ' ' ' ' — 3 



0.39 r 




0.34 E , , , I , , , I , , , I , , , \ , , , 3 

2 4 6 8 10 

smoothing scole R [Mpc/h] 



Fig. 8. — Large-scale asymptotic value of Gnik) as a function 
of the smoothing scale R . T he va lue for R = 1 [Mpc/h], which is 
used for generating Figs. [TO] to [14] is 0.3718. 



with Eqs. (|34p . are the second main results 

of this paper. In the next sections, we shall present the 
detailed assessment of the new terms we have derived in 
this section. 

3.2.2. Shape Dependence 

Let us consider the shape dependence. First of all, the 
last term of Eq. ((32|) has the same shape dependence as 
B^^^, as Qnik) is almost constant on large scale. Thus, 
it peaks at the squeezed configurations as Bj^^ does. 
How about the shape dependence of the other terms? 

All terms in Eqs. ^ have P0(fci) outside 

of the integral, and Eq. ([55)1 contains l/A4R{k) oc fc^^, 
which suggests that all of B^^^, B^F , B]^^\ and 

peak at the squeezed configurations. For sufficiently 
large scales in which T(k) « 1, and for a scale- invariant 
spectrum {P^{k) oc 1/fc ), we may write down Eqs. (|33[) . 
([59]), ((3T|). and jSO]), as 
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i?;f(fci,fc2,fc3)oc('Ji + ^ 

fcl fc2 



5 ^ fcl ■ fc2 f^^^ 
7 2kik2 \k2 fcl 



2 / fcl • fc2 



7 V fc'ifc2 



(2 cyclic), 



|fci - 
|fc2 + g| 



|fc2 + g|^ 
|fci - q\ 



Wnilki - q\)WBiq)T{\ki - q\)T{\k2 + q\. 



xi^2 ("fc25 fc2 + ?) + (5 permutation), 



(40) 



(41) 



Bf{k^,k2,k3)^ 



i3;,«Jfcl,fc2,fc3)cX 



fc2fc3 



2(fc| + fcf) /■ d^q 



|fcl - g|'r(q)W'fl,(g)T(|fci - q\)WR{\h - q|) 



6 



fc2fc3 

+2 



d^g I fcl - g| 
(fcf + fc|) 



1 



1 



fc2 + g|3 

d^g 1 
q |fci - q| 



T{q)WR{q)T{\ki - q\)WR{\ki - q\) 



fca + ql' 

Tiq)WR{q)T{\h - q\)WRi\h - ql] 



fcl - qmq)WR{q)T{\k, - q\)WR{\h - q\'. 



(2 cyclic) 



+ (2 cyclic)(,42) 



(43) 



respectively. 



For a given fci, all of these terms have the largest 
signals when fca is small, i.e., the squeezed configura- 
tions. Note that we do not use the exact scale-invariant 
spectrum for the numerical calculation, but use the 
WMAP 5-year best-fittin g value reported in Table 1 
("WMAP+BAO+SN") of iKomatsu etall (l200l . 
The top-left and bottom-left panels of Fig. [4] show 
and B^^^_^ as a function of fc'2/fci and fca/fci, re- 
spectively, for fcl = 0.01 h Mpc^^. The top- right and 
bottom-right panels of Fig. [4] show the same quanti- 



ties for fcl ^ 0.05 h Mpc" 



We also show i?"*^ 



and 



in the top-left and bottom-left panels of Fig. [5] for 
fcl = 0.01 h Mpc~^, and top-right and bottom-right for 



for fcl 
B^Gi 



In all cases we find that B^^ , 
't ) ^OK,, ^^^'^ peak at the squeezed configura- 

tions, as expected from the above argument. 
We find that the shape dependence of B"^^^ and that 



= 0.05 h Mpc 
and 



of Sy^i, 

/NL ' 



QuG 

9NL 



are quite similar, whereas that of is 

higher toward the elongated triangles, and that of S'iF 

is more sharply peaked at the squeezed configuration. 

We can understand this behavior analytically as fol- 
lows. In order to simplify the analysis, we consider a 
scale- invariant curvature power spectrum , = Pcj,o/k^, 
on large scales where Eqs. dH]), gO]), (glf, and ((43]) 
are valid. On such a large scale, MR{k) can be ap- 
proximated as MR^k) ~ 2fc^/(3iJQrim) = Mok'^, where 
Mo =i 2.16 X 10"^ (0.277/rJ™) [Mpc//i]2 is a constant. We 
focus on the squeezed triangle, fci — fc2 — ak^ = k 
{a ;2> 1), where the signals of all the primordial non- 
Gaussianity terms are maximized. The triangles in this 
configuration lie on the upper side of the triangular re- 
gion of (fca/fci, fc2/fci) plane in Fig. [4] and Fig. [5l and the 
triangle approaches the exact squeezed limit as a — *■ oo. 
With this parametrization, we compare the dominant 
contributions of each of these primordial non-Gaussianity 
terms. 

First, we shall analyze the terms proportional to /nl^ 
B"'^°, Blf, and B"^^. The largest contribution to B"^° 

/NL ' ™ ' /NL " /NL 

in the squeezed configurations occurs when k^ is in the 



denominator: 



oriGO _ n 1^3 pZ 



:4aAlgP|o 

:>nG 



fci 



_fc^ 

fc2fc3 
fc2fc3 fc3fcl 



fc| 



fcafci fcifcj 



(44) 



To compute B^^ , which contains Fj"'', we note that, 
in the squeezed limit, the angular cosines between two 
wave vectors are fci ■ fc2/(fcifc2) = — 1 + l/a^ ~ —1 and 
fc2 • fc3/(fc2fc3) = fcl • ks/ikiks) = -l/(2a). Wc thus find 



b: 



1 



a 



1 



1 



a 



1 



26 

y 



3 d2 



14a2_ 
(45) 



The detailed analysis for B^^^^ is more complicated, as 
Eq. ([55)1 involves a non-trivial integration. We simplify 
the situation by only analyzing the dominant term, which 
can be written as 



B 



■nGl , 
/nl 



'''■H(fcl,fc2) + ^iH(fc2,fcl) 



SaMlP' 



h fcs 

[H(fcl,fc2)+H(fc2,fcl)], (46) 

where 7i(fci,fc2) is the integration that appears in 
Eq. including l/(2cr|j) pre-factor. Note that this in- 
tegration depends only on the magnitudes of two vectors 
and the angle between them; thus, for the squeezed con- 
figuration we are interested in here, 7Y(fc2,fci) depends 
only weakly on a - they depend on a only through the 
inner product of fci ■ fc2 = fc^(— 1 -I- 1/a^). 

Second, we analyze ■ We find that the first cyclic 
terms in Eq. (|30[) are small in the squeezed limit, and the 

is given by 

X(fci) , X(fc2) , X(fc3)" 



dominant contribution to Pg^ 



P"«, = 12M^P 



00 



■■i^MlP^o 



fc2fc3 

2:(fci) 



fcsfci 

2:(fc2) 



fclfc2 



fcofc: 



2^-3 



fcgfcl 



(47) 
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Fig. 9. — Shape of the integration that appears in the dominant 
term of B"f and , Eq. We use four different smoothing 

scales: R=l ,2, 5, and 10 Mpc/h. 



where we have defined 



1 



^Mniq)Mn{\k 



q\)PM- (48) 



We find that T(k) ~ 0.5 and is almost independent of 
k on large scales (e.g., k < 0.03 h Mpc~^ for R 



1.0 Mpc/h; see Fig. IH]). Therefore, by writing X{k) 
we obtain 



B 



nG 



2o, 



(49) 



These results show that all the terms we have analyzed 
analytically so far, B^j^^ , B^f , B]^^, and B'^^^, have the 
same shape (i.e., a) dependence in the squeezed configu- 
rations: they both increase linearly as a increases. This 
explains the shape dependence computed from the full 
numerical calculations presented in Fig. [3] and the top 
panels of Fig. [5] 



Finally, we analyze B 



nG 
f2 ■ 



We find that the second 



cyclic terms in Eq. (|3ip are small in the squeezed config- 
urations. The dominant terms are: 

B^l^SMlPlo 



fcg 



kf 



kikskl 
kf 



kl 



kik2kn 



1,3 



:>nG 



(50) 



Therefore, B^2 increases more sharply as it approaches 
the squeezed limit, oc o;'^. 

This sharp increase of i?'/ F relative to the other terms, 

and that there are many new terms that are of the same 
order of magnitude as im ply that the formula de- 

rived by ISefusatti fc Komatsul (|2007t ) may not be valid 
in the squeezed configuration, where S^F may domi- 

nate over Bj^^ . This is particularly important because 
it is the squeezed configuration that gives the largest sig- 
nal from the primordial non-Gaussianity. We shall study 
this point in more detail in the next section. 



A careful inspection of Eq. p2|) shows that the second 
term within the square bracket diverges when fc2 + 9 = 
or fcs -|- q = 0. This is due to the fact that Pcf,{k) oc 
l/fc^~"= and thus ^^(A:) diverges as fc ^ for Hs < 4. To 
avoid the divergence we set P(f,{k) ~ at k < kmim and 
use kmin = 10~^ h Mpc~^. Fortunately the divergence 
is mild and the results on the squeezed configurations, 
for which i??F gives the most important contribution, 

are insensitive to kmin- changing kmin = 10~^ h Mpc""'^ 
to kmin = 10~^ h Mpc~^ rcsults in negligible changes in 
the squeezed configurations. 

On the other hand, the folded and equilateral configu- 
rations are more sensitive to kmim and we find that the 
difference between kmin = 10~^ h Mpc""'^ and kmin = 
10^^ h Mpc~^ is scale-dependent: at ki = 0.01 h Mpc~^ 
the differences are negligible for all shapes, whereas the 
differences reach ^ 40% at k ^ 1 h Mpc~^. (Note that 
the difference in the squeezed configuration reaches 1% at 
k ^ 1 h Mpc~^, being totally negligible on larger scales.) 
While this divergence does not have much observational 
consequences (because the signals in the other configu- 
rations at A: > 0.01 h Mpc~ would be dominated by the 
other non-linear effects: B^, non- linear bias and terms 
proportional to /nl, as we show in the next section), 
there may be a better treatment of the divergence than 
setting Piji{k) — at k < kmin- 

3.3. Scale Dependence 

How important are the primordial non-Gaussianity 
terms, B?*^", B^^i, B^f , S"F , and B""^ , relative to 

' JNL ' JNL ' ™ ' /^j^ ' 3NL ' 

B^ and the non-linear bias term? Which one is the most 
dominant of the primordial terms, terms proportional to 
/nl, , or ? How about the scale-dependence? 



How about the shape dependence? 
We collect all the terms proportional to /nl, and call 

it 



Tjtot 7^nC?0 

JD f =L JD £ 
/NL jNL 



62 



B 



nGl 



(51) 



Throughout this section, we use the standard value of 
62/61 = 3(127r)2/3/20 ~ 1.68 from a spherical collapse 
model. 

Figure [10] shows the scale and shape dependence of 
each term in Eq. ((5T|) evaluated at z = 0. For all configu- 
rations shown in thi s figure, the primordial non-G aussian 
term calculated in ISefusatti fc Komatsul (j2007l ) is the 
smallest among four /nl terms, which means that the 
non-Gaussian signal on large scales is much bigger than 
recognized before. 

For the squeezed triangle, all of the terms in Eq. (|5ip 
depend on ki in a similar way. We find their ratios by 
comparing Eqs. (O, (gS]), and (gH): 



B 



nGO 
/nl 



B 



nG 



B 



nGl 
/nl 



1:^ 

28 



2.96. 



(52) 



Note that we have used the numerical value of 
TLiki^kq) ~ 0.741 for a = 100, and this value slightly 
increases when a decreases^. Therefore, for the squeezed 

® On large scales, k < 0.01 h/Mpc, the numerical ratio 
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Fig. 10. — Scale and shape dependence of the galaxy bispectrum terms that are linearly proportional to /nLi as a function of fci. 
Except for the bottom-right panel, we use R = 1 Mpc/h. (Top Left) The squeezed triangles with fci = ^2 = lOOfca, (Top Right) 
the elongated triangles with ki = k2 + fca and k2 = 3^3, (Middle Left) the folded triangles with ki = 2k2 = 2^3, (Middle Right) 
the isosceles triangles with 3fci = 4^2 = 4fe3, and (Bottom Left) the equilateral triangles with fci = A;2 = ^3. The thick dot-dashed, 
dashed, solid, and dotted lines show the contributions from the primordial non- Gaussianity: the 3"^^° (Eq. b2/HB^ (Eq. |[33t). 



4(62/^1) [6ij(fcl) + Gfi(fc2) + eH(fc3)] -B^j^j^^ {GRik) defined in Eq. and b2/biB'^°^ (Eq. ^) terms, respectively. The thin dotted 

and dashed lines show the non-linear effects: (Eq. ^21) ) and the non-linear bias (the second term in Eq. respectively. We use 

the standard value of 62/^1 = &c — 1.686 from spherical collapse model. (Bottom Right) Dependence of the squeezed bispectrum on the 
smoothing scale, R, showing that the dependence is negligible for k\ <C 1/R. 



triangle, we find a simple and illuminating result: 

fc2, fcs) - ISSy^f (fci, fc2, fca). (53) 
This is an important result, showing that the statistical 
error on /nl from the galaxy bispectrum will be smaller 
by at least a factor of 15, comp ared to what was obtained 
in lSefusatti fc Komatsul ()200'it ). 

Figure [TT] and the top panels of Figure [12] show vari- 
ous bispectrum terms in various triangle configurations 
(see Fig. [1] for the visual representations of the triangles), 
evaluated at z = 0. As an example we use the follow- 
ing bias and non-Gaussianity parameters: 62/^1 = 0.5, 
/nl = 40, and (7nl = lO**. The value of the linear bias, 
61, is irrelevant here as it does not change the relative 
importance of terms in Eq. (|14p . and thus we show the 
bispectrum terms divided by h\. 

The message is quite simple: it is the squeezed configu- 
ration that provides the best window into the primordial 

^/i?L /^/i?L constant, and is equal to 3.15, 3.06, 3.00, and 2.98 
for a = 10, 20, 50, and 100, respectively. 



non-Gaussianity. The other non-linear effects become 
more and more dominant as we move from the squeezed 
to the equilateral, i.e., (a) to (c) in Fig.[TJ Even with this 
generous amount of non-Gaussian signals, /nl = 40 and 
ffNL = 10'*, only /nl term can be visible in the isosceles 
and equilateral configurations on large scales. 

For the the non-squeezed configurations, the /^^ ^"^^ 
gNL terms with the above chosen parameters are com- 
parable and the /nl term is order of magnitude greater 
than the /^l ^^'^ ffNL terms; however, the /^^ term is 
the most dominant of all on large scales in the squeezed 
configuration [a > 10). 

We can understand these results analytically by com- 
paring Eqs. dSSl), (gH), (gni), and dSO]). For the squeezed 
triangles with fci = fc2 = ak^ (a 3> 1) and a scale- 
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Fig. 11. — Scale and shape dependence of various bispectrum terms, B{ki, k2, k^), as a function of ki. For the figure except for 
the bottom right, we use R = I Mpc/h. (Top Left) The squeezed triangles with fci = ^2 = lOOfcs, (Top Right) the elongated triangles 
with fci = fc2 + fcs and k2 = Sfca, (Middle Left) the folded triangles with ki = 2k2 = 2^3, (Middle Right) the isosceles triangles with 
3fei = 4^2 = 4fc3, and (Bottom Left) the equilateral triangles with ki = k2 = k:j. The thick dot-dashed, triple-dot-dashed, and solid 
lines show the contributions from the primordial non-Gaussianity: the ft^i^B^^^ (Eq. I ISll l). fe/^iSNLSgj^L C^l- Gnj): ^^^d ''a/fei/^L^j? 
(Eq. (|3TJ) terms, respectively. The thin dotted and dashed lines show the non-linear effects: (Eq. I I21I I') and the non-linear bias (the 
second term in Eq. JOj), respectively. We use the standard value of 62/^1 = 5c — 1.686 from spherical collapse model. (Bottom Right) 
Dependence of the squeezed bispectrum on the smoothing scale, R, showing that the dependence is negligible for ki <^ 1/ R. 



invariant spectrum, (y: k ^, wc find 



/nlM°1 15 M^e 



^ 0.0240 



100 



40 



a ) /nl VO.Ol h Mpc" 



- (54) 



gNL-Bgi^ 3 ffNL 

f2_ ^« JnL 



f2 n,iG 



-0.000938 



-25 6^— 

40 5NL VO.Ol h Mpc"\ 



(56) 



These estimates confirm that B f2 dominates over Bf.,, 

JNL J^^^ 

and Bgj^^ in the squeezed configurations on large scales, 
k < 0.05 h Mpc"^ for a ^ 100, and /fc < 0.03 h Mpc"^ 
for a = 50. For a = 10, /^l term dominates only on the 
extremely large scales: k < 0.006 h Mpc~^. 
Note that for a given configuration (for a given a) 



B"^ /B 

/NL' 



nG 



P and B"'^ /B"G oc fc^ while B"'^ /B"? 

jNL' SNL SNL' /j^l 

is independent of k, which is consistent with what we 
show in Fig. [11] on k < 0.1 h Mpc"^ 

In summary, the most unexpected and important re- 
sults of our study are as follows. 

• The ter ms that are linearly p roportional to /nl, de- 
rived in lSefusatti fc Komats u (2007), receive addi- 
tional contributions, and are enhanced by a factor 
of 15 for the squeezed triangles (see Eq. ((5^ ). 

• The (oi' ■''nl) term actually dominates over the 
/nl term by a large factor for the squeezed triangles 
(see the top- left panel of Fig. [TT|) . 
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111! (Middle) 2 = 3 and b2/h\ = 1.5. The non-Gaussianity parameters, /nl = 40 and = 10*, are the same as in Fig. 1111 
2 = 3 and b2/b\ = 1.5. The non-Gaussianity parameters, /nl = 4 and jnl = 100. 



This suggests that the galaxy bispectrum is more 
sensitive to /nl t han previously recognized by 
iSefusatti fc Komatsul (|2007f ). greatly enhancing our abil- 
ity to detect the primordial non-Gaussianity of local 
type. On very large scales, fci ^ 0.01 h Mpc~^, even 
the _gNL term (with ^nl = lO'*) dominates over the /nl 
term, giving us a hope that perhaps we can obtain a 
meaningful limit on this term using the galaxy bispec- 
trum. 

3.4. Redshift Dependence 



The quantities we have calculated so far are evaluated 
at the present epoch, z = 0. At higher redshift, each 
quantity needs to be scaled with some powers of the 
linear growth factor D(z), which is normalized to 1 at 
the present epoch. We find Pr oc D'^{z), B^^ cx D'^{z), 
B'f° cx D^{z), B^'^ oc D^{z), B"}^^ oc D^{z), S'/F oc 

D^(z), and B^^^ oc Bi^{z). Therefore, the final result for 
the halo bispectrum from the local type non-Gaussianity 
is 



Bg(hrMM.z) = h\[z)D^{z) 



Bg{ki,k2, fe) + ^ {PR{ki)PR{k2) + (2 cyclic)} + J^B]^°{ki,k2, fe) 



h{z) 



D{zy 



Mf) / (B^fiki, fc2, fe) + 4 (GRiki) + gR{k2) + GRiks)) B]g°{ki,k2,k3) + B]^^\ki,k2, ks)) 



bi{z) \D{z) 



(57) 



where B^, Pr, B^j^^ , Bl^ , -^/f ■ B^^^ arc From equation ^Ji it is clear that the contributions 

evaluated at z = '^^ from non-Gaussian initial conditions become more and 
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more important as wc go to higher redshifts. The new 
terms that we have derived in this paper, the B"p and 

-^ffNx terms, are even more important than the term de- 
rived bv lSefusatti & Komatsd (I200I . 3"}^°, This prop- 
erty makes high-redshift galaxy surveys particularly a 
powerful probe of primordial non-Gaussianity. 

Fig. [T3] and the middle panel of Fig. [T^ show the bis- 
pcctrum terms at z = 3. Note that wc use a larger value 
for the non-lin ear bias, 62/61 = 1.5, in accordance with 
a halo model (jCoorav fc Shethll200^ . At this redshift, 
with /nl — 40 and ^nl = 10^, the ^nl and /^^ terms 
dominate over the non-linear effects also in the elongated, 
folded and isosceles configurations at fc < 0.01 h Mpc^^, 
as well as in the squeezed ones. The /nl terms dominate 
over the non-linear effects on even smaller scales, and the 
importance of the /^^ a-nd 5nl terms relative to the /nl 
term is greater, as expected from their dependence on 
D{z). 

4. DISCUSSION AND CONCLUSIONS 

Let us come back to this question, "can we still use 
Sefusatti & Komatsu's equation, Eq. ([T]), with bi re- 
placed by the scale-dependent bias, Eq. (jl])?" The an- 
swer is clearly no: the primordial non-Gaussianity gives 
the largest signal in the squeezed limit, whereas the non- 
linear gravitational evolution and non-linear bias give the 
minimal signals in the same limit. This means that these 
effects are physically totally distinct, and thus a mere 
scale-dependent rescaling of one effect docs not give an- 
other. Therefore, replacing bi in Eq. ([T]) with the scale- 
dependent bias in Eq. ^ results in an incorrect predic- 
tion. For example, even though we have a term similar 
to that of the scale-dependent bias, , in our final 
expression of the galaxy bispectrum for the local-type 
primordial non-Gaussianity, there are many more terms 
that do not look like the scale-dependent bias that ap- 
pears in the galaxy power spectrum. Furthermore. 
is by no means the most dominant term. 

In this paper, we have derived a general expression for 
the bispcctrum of density peaks in the presence of pri- 
mordial non-Gaussianity (Eq. using the MLB for- 
mula as well as using the local bias ansatz. This result 
is general as long as we consider the bispectrum of high 
density peaks, i.e., v = Sc/cfr 3> 1, which is equivalent 
to highly biased galaxy populations, 61 ^ 1, on large 
scales in which the n-point correlation functions are much 
smaller than unity. (This condition was necessary for us 
to Taylor expand the exponential in Eq. (O.) 

Wc have applied our formula to the local form of pri- 
mordial non-Gaussianity in Bardeen's curvature pertur- 
bations, ^ = (j) + f^Lff)'^ + ffNL^'^, and found new terms 
that are proportional to /nl, /mt and gNL, which were 
absent in the formula derived by ISefusatti fc Komatsd 
()2007| ). We have examined the shape and scale depen- 
dence of these new terms as well as those of the known 
terms, and found that the primordial non-Gaussianity 
contributions yield the largest signals in the squeezed tri- 
angle configurations, where the non-linear gravitational 
evolution and non-linear bias yield the minimal signals. 
This is a good news: this property enables us to distin- 
guish the primordial and non-primordial effects easily. 

The effects of primordial non-Gaussianity on the 
galaxy bispectrum arc more important in a high red- 



shift universe, and thus high-redshift galaxy surveys are 
particularly a potent probe of the physics of inflation via 
measurements of primordial non-Gaussianity. 

The most significant conclusion of this paper is that, 
in the squeezed configurations, the /^^ term actually 
dominates over the /nl term by a large factor, and, 
on large scales, newly derived /nl term dominates over 
the non linear terms for all configurations. Because of 
this, the galaxy bispectrum should be more sensitive to 
/nl than previously recognized: we have found a fac- 
tor of ~ 15 enhancement f or the /nl term studied in 
ISefusatti fc Komatsd ()2007t ). In addition it is also sen- 
sitive to a new term, pnl- Figure [Ml and the bottom 
panel of Fig. [T^ shows the bispectrum at z = 3 with 
much reduced primordial non-Gaussianity parameters, 
/nl = 4 and pnl = 100. In the squeezed configurations 
{a = 100), the /^^ term is still well above the usual terms 
from Gaussian initial conditions at fc < 0.1 h Mpc~^, the 
/nl term is above at fc < 0.4 Mpc^^, and the (7nl term 
is above at fc < 0.01 h Mpc~^. Even with the milder 
squeezed limit for a — 10, the /nl term still is above the 
Gaussian term at fc < 0.02 Mpc~^. 

The fact that the /^^ term dominates in the squeezed 
limit is particularly interesting, as it provides us 
with the unique window into the physics of infla- 
tion in the following way. Recently , a number of 
group s (e.g.. IBoubekeur fc LvthI [2009: iHuang fc Shid 
I200a iBvrnes et all l2006t iBuchbinder et al have 
shown that the primordial trispectrum can in general be 
written as 



T^{ki,k2, ks, ki) 



25 



= 6gNL [P4>{ki)P4>{k2)P4,{k3) + (3 cyclic)] + — tnl 

io 

X [P4,{ki)P^{k2) {P^ikis) + P^{ki4)} + (11 cycliO§§) 

instead of Eq. (fTT]) . Different models of the early uni- 
verse predict different relations between tnl and /nl- 
Therefore, separately detecting the tnl (i-e., /nl) and 
/nl terms can be a powerful tool for constraining the 
model of the early universe. 

How well one can constrain these parameters with the 
planned future high-redshift galaxy surveys will be pre- 
sented in a forthcoming paper. 
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APPENDIX 
A. INTEGRATION OF r^"^ 

In the standard perturbation theory, the four-point correlator contained in the definition of T^^^ (see Eq. ([^5])') is 
given by^ 

(<5(i)(fci)<5(i)(fc2)<5(i)(fc3)<5(2)(fc4)) = J ^F^'\q,h - q){6<^'Hki)S^'\k2)d^'Hh)S^'\h - q)6^'Hq)). (Al) 

For non-Gaussian density fields, the leading order of Eq. (jAip contains the ensemble average of products of six Gaussian 
variables, </>, which gives products of three power spectra, P^. We find 

d3q 



(27r)3 2 

W) 



F^'\q, ki - (fc2)<5(^) ih)6^'^ {h - q)S'^'Hq)) 

q){5^'Hki)5^'\k,)){6'-'\h)d^'\h - q)S^'Hq)) + (cyclic 123) 



F^'\q,h 



3^2 



^F^'\q,h - q){S^'Hk,)6^'Hq)){S^'\k2)S^'Hh)5^'\k, - q)) + (cyclic 123) 

--{27Tf 2/NLP,n(fci)X(fc3) J d\M{q)M{\k^- q\)P^{q){P^i\ki- q\) + 2P^{k3)}F^'\q,k4- q)6''{h2) 
+Af^LM{k2)M{ks)M{ki4)P^{k,)F.^'\~kuku) 



X {P4k2)P4k3) + P^{k2)P^{ku) -f P^{k3)P^{ku)} + (cyclic 123) 

® ISefusattH I I2009I) also derived and studied this term independently. 



5''{ki2Zi)- 



(A2) 
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Therefore, T^^^ is given by 



Tk'^\ki,k2,k3,ki) 

--WR{kl)WRik2)WRik3)WRik4) 



2hLPmiki)Mik3) J d-'qM{q)M{\ki - q\)PM {^0(1^4 - q\) + 2^0(^3)} F!f\q, h - q)S''{ki2) 

+4hl^M{k2)M{k3)M{ku)Prn{kl)F^'\^h,ku) 

X {P4k2)P^{k3) + P0(fc2)P0(fci4) + n (^-3)^0(^14)} + (cyclic 123) 
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where A4{k) = A4u{k)/Wii{k), kij = ki + kj, and (cyclic 123) den otes that the cychc changes among {ki, ^2, fca). We 
calculate the sum of {1112} terms in Eq. ([27|) by integrating Eq. (jA3|: 



J (27r)3 

8hLWRik2)WR{k3)M{ki)M{k2)M{k3) {P^ik2)P4k3) + P^(fc2)P0(fcl) + ^^(fcs)^^^!)} 

X y ^#fl(|fci-g|)iyflX'7)^m(g)^^i'^(-g, fci) 

-4.hLWRik2)WRik3)Miki) [Pmik2)F^'\k2, h) + P„,{k3)F^'\k3, h) 

-^WRilh - q\)WR{q)M{q)M{\h - q\) {P^mkl ~ q\) + 2P^{q)P^(kx)} 



{2n) 

+ &hLWRik2)WRik3)M{k3)P„,{k2) J ^WRilh - q\)WRiq)Mi\ki - q\)M{\k2 + q\)F!f\-k2, fca + q) 
X {P^{k:i)P^[\ki - q\) + P^{k^)P^{\k2 + q\) + P^{\ki - q|)P^(|fc2 + q\)} + (fca ^ h) 

+ 8/NLM^«(fc2)T^i?,(fc3)X(fc3) J -0j3^Ri\h - q\)WRiq)M{\h - q\)Mi\k2 + g| )P™ (g)^^-' ^ 9, ^ + q) 
X {P^dfci - q\)P^ik3) + P^dfci - q\)P4\k2 + q\) + P^ik3)Pc^{\k2 + q\)} + (fc2 ^ h) 

+ 8/NL (l^fl(^-2)VKfl(A:3))' Prn(^-3)A^(fc2) 

X J ^Miq)Mi\k2 - q\)PM {^4|fc2 - q\) + 2P^(fc2)} P^'^ - g) + (fe ^3). (A4) 

B. SUMMARY OF EQUATIONS 

As various terms contributing to the galaxy bispectrum are scattered over various places in the paper, we collect 
them together in this Appendix. For galaxies of size R (or mass M ~ (47r/3)P^pm, where p,„ is the cosmic mean 
matter density), the galaxy bispectrum at redshift z is given by 



Bg{ki,k2M.z) = hl{z)D\z) 



Bgik,, fc2, fcg) + ^ {Pfl(fci)Pfl(fc2) + (2 cyclic)} + J^B'}^^^{k,,k2, ks) 



(^"""(^i' ^3) + 4 {GRik,) + gR{k2) + GRik,)) P7«0(fci, fe, fca) + P;,T(fci, fc2, k,)) 



^p;£(fci, fc2, ^3) + ^f'^olik.M. k,) 



(Bl) 



where bi{z) and b2{z) are the linear and non-linear bias parameters, respectively. As we mentioned in §2, b2{z)/bi{z) 
would be equal to 6c within the context of the MLB formalism, but the precise value has to be measured from N-body 
simulations. 

Note that the redshift evolution of each term in explicitly given by the powers of the linear growth factor D{z), and 



Primordial non-Gaussianity and galaxy bispcctrum 
various contributions, B^, Pr, B"^ , B'^^°, B^j^^ , BjF , and B^^^, are evaluated at z = with 

Bg{ki,k2, h) = 2F^'\kuk2)WRih)WR{k2)WRih)P^{ki)Prn{k2) + (2 cyclic) 



2^iM^iMA,,(fc3) + (2cyelic) 



^MR{ki)MRik2 

Blf{ki,k2. k3)=AWR{ki)WR{k2)WRik3) 

1 



^niki) J^R{k2) 



Pm{ki)P,n{k2)F^'\ki,k2) + (2 cyclic) 



MRik,) MR{k2) 

8WRik2)WRik3)Pm{k2)M{k3)P^{k3) J vS^W^fld^i " q\)WR{q)M{\ki - q\)Mi\k2 + q\) 



[P4,i\k2 + q\) + P4\ki - q\)] F!f\-k2, k2 + q) + (5 permutation) 
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(B2) 
(B3) 

(B4) 
(B5) 



iiMR{k2)MR{k-i)P^{ki) [P^{k2) + P^lfcs)] 



(2^) 



^A^fl(g)A1ii(|A;i - q\)P^{q) + (2 cyclic 



Rik2)MRik3)P4k2)P4k3) J ^MRiq)MRi\h - q\ 



B^gik,,k2,k,) 



X [P^dfcs + q\) + P^dfe + g|)] + (2 cyclic) 
1 



(B6) 



6MRik2)MR{k3) [P4fc2) +P0(fc3)] y j^MR{q)MRi\ki - q\)P^{q)P^{\ki ~ q\) + (2 cyclic) 



- 12A^fl(fc2)Xii(fc3)P4fc2)P^(fc3; 



(2^ 



- 91)^(9) + (2 cyclic) 



(B7) 



Note that we show only dominant terms for B'^i'^'^ and P^F on large scales. One can find the exact definitions in 
Eq. dH]) and Eq. ([M]). Finally, J^/f(fc) and QR{k) are defined as follows. 



TR{k)^^^ J ^P^{q)MR{q)MR{\k-q\ 



PA\k-q[. 



^-(^)"2^i(2^ ^^^^(^^ ^"('^^^^ 



(B8) 
(B9) 
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